Astronomy & Astrophysics manuscript no. Das 


February 2, 2008 


(DOI: will be inserted by hand later) 





Computation of outflow rates from accretion disks around black 

holes 

Santabrata Das, Indranil Chattopadhyay, A. Nandi, and Sandip K. Chakrabarti* 
S. N. Bose National Centre for Basic Sciences, Salt Lake, Kolkata 700098, India 

Received ISth May 2001 / Accepted SOth August 2001; Astronomy and Astrophysics, 379, 683, 2001 

Abstract. We self-consistently estimate the outflow rate from the accretion rates of an accretion disk around 
a black hole in which both the Keplerian and the sub-Keplerian matter flows simultaneously. While Keplerian 
matter supplies soft-photons, hot sub-Keplerian matter supplies thermal electrons. The temperature of the hot 
electrons is decided by the degree of inverse Comptonization of the soft photons. If we consider only thermally- 
driven flows from the centrifugal pressure-supported boundary layer around a black hole, we find that when the 
thermal electrons are cooled down, either because of the absence of the boundary layer (low compression ratio), or 
when the surface of the boundary layer is formed very far away, the outflow rate is negligible. For an intermediate 
size of this boundary layer the outflow rate is maximal. Since the temperature of the thermal electrons also decides 
the spectral state of a black hole, we predict that the outfiow rate should be directly related to the spectral state. 
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1. Introduction the viscous Keplerian disk resides in the equatorial plane, 

•"^ r while the weakly viscous sub-Keplerian flow flanks the 
■ Most of the galactic black hole candidates are known to 

5^ ■ ^ ^ ^ Keplerian component both above and below the equatorial 

' undergo spectral state transitions (Tanaka & Lewin, 1995; 
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plane. The two components merge into a single compo- 
nent when the Keplerian disk also become sub-Keplerian. 
It is suggested (Chakrabarti, 1990) that close to a black 
hole, at around 10 — 15 Vg, [vg — 2GMbh I (? is the 
Schwarzschild radius, Mbh and c are the mass of the 
black hole and the velocity of light respectively) the sub- 
Keplerian flow slows down due to the centrifugal barrier 
and becomes hotter. Chakrabarti (1999, hereafter Paper I) 
shows that this centrifugal pressure-supported boundary 
layer (CENBOL for short) region could be responsible for 
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Chakrabarti & Titarchuk, 1995, hereafter CT95; Ebisawa 
et al. 1996). Two common states are the so-called hard 
state and the soft state. In the former, soft-X-ray lu- 
minosity is low and the energy spectral index a ~ 0.5 
{E^ cx i^^") in the 2-lOkeV range. In the latter state, the 
soft-X-ray luminosity is very high, and hard-X-ray inten- 
sity is negligible. There is also a weak power-law hard- 
tail component with an energy spectral slope a ^ 1.5. In 
the two component advective flow (TCAF) model (CT95), 
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jets and computed the ratio of the outflow to the inflow ordinate in which the black hole is treated at the cen- 
rate assuming a simple conical accretion disk model. tre. We use geometric units in which all the length, time 
In the present paper, we compute the absolute value of and velocity scales are measured in units of 2GMbh/c^-, 
the outflow rate as a function of the rates of the two in- 2GMbh/c^ and c respectively. In future, we use r to de- 
flow components, Keplerian and sub-Keplerian. This we note non-dimensional distance, d and a to denote the non- 
do analytically following the recently developed proce- dimensional radial velocity and adiabatic speed of sound 
dure of obtaining shock locations (Das, Chattopadhyay respectively. In accretion or outflow, we assume that the 
and Chakrabarti, 2001). By dynamically mixing these two viscous stress is negligible so that matter moves with a 
components using solutions of the viscous transonic flows constant specific angular momentum. Indeed, even if vis- 
we obtain the specific energy and angular momentum of cosity is not negligible, the transport of angular momen- 
the sub-Keplerian region. We use these pair of parameters tum is slow compared to the infaU timescale. Hence, mat- 
to locate shocks in the flow, compute the compression ratio ter can have almost constant specific angular momentum, 
and from this, the outflow rate. We note that as Keplerian In this case, the radial momentum equation for a 
matter is increased in the mixture, the shock compression non-dissipative flow in vertical equilibrium is given by 
ratio goes down, and the outflow rate decreases. This is (Chakrabarti, 1989), 
also the case even from a radiative transfer point of view (I'd 1 dP 1 

'^ir + --i — ta2=o- (1) 

- when the Keplerian rate is high, the CENBOL region is ^ ~ ^1 

completely cooled and the shock compression ratio I. Integrating this, we obtain the conserved specific energy 

Hence in the soft state, which is due to increase of the of the flow, 

Keplerian rate, outflow should be negligible. 1 „ n }? 1 

£, = -i} +na + _-___, (2) 

In the next Section, we present the governing equa- ^ ' 

tions to compute the outflow rates using a purely analyt- ^here n is the polytropic index of the inflow and A is 

ical method. We compute resuhs for both the isothermal t^^c specific angular momentum. In Eq. (1), P and p are 

and adiabatic outflows. In §3, we present our resuhs for a thermal pressure and density respectively, v is the infah 

single component sub-Keplerian flow. We also produce ex- velocity and a = ^/{jP/p) is the adiabatic sound speed. 

amples of realistic disks with Keplerian and sub-Keplerian The mass flux conservation equation in a flow which is 

components and obtain outflow rates as functions of the ^ vertical equilibrium is given by, 

inflow parameters. In S4, we discuss our results and draw ,v n i 

^ Min = A-Kfy&rh{r) = Qinps'&srl, (3) 

conclusions. 

where Qin{— ^J^^^T^CLsfV^) is the solid angle subtended 
by the inflow at the CENBOL boundary. Subscripts s de- 

2. Model Equations 

notes the quantities at shock (CENBOL boundary) and 

We consider matter accreting on the equatorial plane of a h{r) = ^J^ar^/^{r — 1) is the half-thickness of the disk in 

Schwarzschild black hole. Spacetime around the black hole vertical equilibrium at a radial distance r. 
is described by the Paczynski-Wiita pseudo-Newtonian A sub-Keplerian flow with a positive energy will pass 

potential (j) = r-'iGM^Jn/c'^ (Paczynski & Wiita, 1980) through the outer sonic point and depending on whether 

where Mbh is the mass of the black hole and G, c are the the Rankine-Hugoniot condition is satisfied or not, a 

gravitational constant and velocity of light respectively, standing shock may form (Chakrabarti, 1990; Chakrabarti 

Here, r is the radial distance from the origin of the co- 1996). If a standing shock forms, then the post-shock 
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flow would become hotter and would emit hard X-ray ra- sign was written using the mass flux conservation equation 

diation. This CENBOL region behaves similarly to the given above (Eq. 3). 

boundary of a normal star; it would be expected to drive At the shock, the total pressure (thermal and ram pres- 

outflows. Using Eqs. (2) and (3), it is easy to obtain shock sure) is balanced: 

locations (i.e., outer surface of the CENBOL) analytically. ^ i:_{r,)^^_{r,) = W+{rs) + E+(r,)i?^(r,), (6) 

Briefly, the procedure to obtain shocks involves the follow- 
where W is the pressure of the gas integrated vertically. 

mg steps: 

, , ^ . . r .„ ^ ^ , We assume that in the pre-shock region, the thermal 

(a) i'or a given pair of specihc energy Cy and angufar mo- 

. . . ^ , . pressure is small in comparison to the ram pressure, 

mentum A, one obtains a quartic equation for the sonic 



point and solves it for the three sonic points located out- W+{rs) = — — — 5]_(rs)i?_(rs). (7) 

K 

Side the horizon. Two of them are saddle type or 'X' type r^^^ isothermal sound speed in the post-shock region 
sonic points and one is a centre type or 'O' type sonic -g obtained from- 

r2-A2(r,-l) 



point. = ^ = = - 

(b) From the inner and the outer 'X' type points, Mach * /o 



1) 



(8) 



numbers are expressed as polynomials of radial distance where, fo = -^rj- 

r. These Mach number expressions satisfy constraints that Up to the sonic point matter moves slowly and the 

they must have appropriate values at the sonic points. density is higher. Since the outflow would take place in a 

(c) In addition, it is enforced that the Mach number in- sea of radiation, the momentum deposition is likely to be 
variants at the shock location are also satisfied (r^). effective. With the electron number density Ue oc r~^/^, 

(d) The resulting equation becomes quartic in and the yet photon number density Uj oc r"^, it is easier to de- 
shock locations are obtained from its solution. posit momentum only close to the black hole. In radiation 

Details are discussed in Das ct al. (2001). We consider driven outflows from the stellar surface, it is customary 

only the region of the inflow parameter space {£y, A) that to assume flows to be isothermal until the sonic point. We 

is able to produce standing shocks. first compute outflow rates making this assumption. Later 

In the pre-shock region, matter is cooler and is sub- we drop this assumption and show that the general be- 

Keplerian. Assuming £y ^ (freely falling condition) and haviour remains similar. In addition, we assume that there 

a ^ (cool gas) in presence of angular momentum, matter is very little rotation in the outflow. There is no a priori 

will fall with a velocity, reason to assume this, except that there is no observa- 

■ 1/2 tional support of rotation in the jet and it is possible that 
(4) 

r — 1 J due to radiative viscosity most of the angular momentum 
Using this, from Eq. (3) the density distribution can be is transported very close to the black hole. Furthermore, 
obtained. it has been observed that the effect of angular momen- 
At the shock r = r^, i.e., the boundary of the tum in the outflow is to bring the sonic points closer to 
CENBOL, the compression ratio is given by, the black hole, especially away from the axis (Sakurai, 
S_i_ h^(r )o^(r ) 'd 1985; Chakrabarti, 1986). The general eSect would pro- 
S_ h-{rs)p-{rs) duce a transverse structure in the jet which we ignore in 
where subscripts "— " and refer, respectively, to quan- the present solution. It was shown (Das & Chakrabarti, 
titles before and after the shock. Here, S is the density of 1999) that in presence of angular motion the conical out- 
matter integrated vertically (S ~ ph) and the second '=' flow is to be replaced by an annular flow confined by the 
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centrifugal barrier and the funnel wall. Generally speak- 
ing, the outflow surface varies as r"*/^. However, the in- 
flow surface area is still proportional to . Because of 
this asymmetry, the problem is no longer tractable ana- 
lytically and is beyond the scope of the present paper. 

2.1. When the outflow is isothermal 

The radial momentum balance equation in the outflow is 
given by 

^dP 1 
and the continuity equation is given by 



ar 



{pdr^) = 0. 



(10) 



(11) 



Eliminating ^ from above two equation we get 
M_N_ 

where N = - 2(r-i)'^ and D = 'd- -f. 

To obtain the sonic point condition, we put A'' = and 
D = and get, 'i?(rc) = Cg, and = ^+^'-^t.^V^+i6g, ^ 
where the subscript c denotes the quantities at the sonic 
point in the outflow. 

Integrating the radial momentum equation, consider- 
ing the sonic point condition, we have, 



Cllnp+ 



(12) 



2(r,-l) 2 ^ ' ^ 2(r,-l)- 
Here, we have ignored the radial velocity in the outflow at 
the boundary of the shock. Using the notations p(rc) = pc 
and p{rs) = p+, we obtain, 



Pc = p+exp[-f\, 

where/=i-2^(^^^fy(^. 
The outflow rate is given by 



Mout = QoutPc'&cTc, 



(13) 



(14) 



where Qout is the solid angle subtended by the outflow. 
Prom Eq.(3) & Eq. (14) we get, 



M„ 



Mi, 



e. 



rUrs 



r2-A2(r,-l) 



-1/2 



RCsrj 
rs{rs - 1) 



expl 
(15) 
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The above relation is very similar to that obtained in 
Paper I when the effects of rotation in the inflow were 
ignored. However, there the ratio Rm was a function of 
R alone. In the present analysis, R is computed self- 
consistently from the specific energy and the specific an- 
gular momentum of the flow: 



R 



1/2 



(16) 



where pre-shock and post-shock Mach numbers M_(f , A) 
and M+(f, A) are computed analytically from Das et al. 
(2001). 

2.2. When the outflow is adiabatic 

At the other extreme, when the energy of the outflow docs 
not change, one can also obtain an analytical expression 
for the outflow rate assuming the >> X^. In this case, 
the entropy density of the flow in the post-shock region 
is the same at the entropy density of the entire outflow 
and the spcciflc energy is also conserved along the out- 
flow. We assume that the turbulence generated at the 
CENBOL has effectively transported angular momentum 
away. Thus, the energy conservation equation gives 



1 
2^s 



2n+l 



1 
2^c 



(17) 



where the left hand side is the energy at the CENBOL 
(r = rg) and the right hand side is at the sonic point 
(r = Tc) of the outflow where Uc = Uc has been used. 
n = is the polytropic constant. In a Bondi (in or 
out)flow, = i/4r-^. At the CENBOL, = jC^, where 
Cg is the isothermal sound speed (Eq. 8). Using these, one 
obtains (assuming >> A^) 

re ^ 2n-3 

r, 4(^-1)' 

and 

In an adiabatic flow with an equation of state P = Kp'^ 
(where K is & constant and a measure of entropy), one 



obtains, assuming, Kc = Kg , 

Ps 



(19) 
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From these relations one obtains the ratio of the outflow 
to the inflow rate as 

0o,/o,3-Rr4r8(^-l) 



Rri 



i?2 



- (20) 



Here, we have used n = 3 for a relativistic flow. The nature 
of this function will be discussed below. 

3. Outflow rates from inflow parameters 

In Eq. (15), we presented the outflow/inflow rate ratio 
as a function of the compression ratio of the flow at the 
shock. The compression ratio is obtained from the specific 
energy and angular momentum using Eq. (16). First, we 
employ analytical means to obtain this for a single compo- 
nent sub-Keplerian disk. Second, we use a two component 
Keplerian/sub-Keplerian disk to actually compute these 
parameters from more fundamental parameters such as 
accretion rates and viscosity. 

3.1. Single component sub-Keplerian Flows 

In Fig. 1, wo plot the analytical solution of the compres- 
sion ratio i? as a function of the flow parameters: spe- 
cific energy £ and the specific angular momentum A. The 
shock strength generally increases when energy decreases 
and the angular momentum increases. This is because for 
low energy, the outer sonic point and the shock form very 
far away and the Mach number jumps from a very large 
number to a very small number. If the angular momen- 
tum is decreased, shock is produced only if the specific 
energy is high, i.e., if the sonic points and the shocks are 
very close to the black hole. Here, flow becomes subsonic 
before its Mach number could be very high. 

Figure 2a shows the principle result of our work when 
only one sub-Keplerian accretion is chosen as the inflow. 
We plot the ratio Rm for a large number of specific angu- 
lar momenta of the flow ranging from 1.57 (innermost) to 
1.83 (outermost) at intervals of dX = 0.02. The curves are 
drawn for aU ranges of specific energy £ for which shocks 
are formed. Along the X-axis the compression ratio R of 
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Energy (S) 

Fig. 1. Variation of the compression ratio of the shocks 
as a function of specific energy £ and angular momentum 
A as obtained from the analytical solution. A varies from 
1.57 (right) to 1.79 (left). Curves are drawn at intervals of 
dX = 0.02. 
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Compression Ratio (R) 

Fig. 2a. Variation of the ratio of outflow to inflow rates R,h 
as a function of compression ratio for various specific angular 
momenta. A = 1.57 (inner most) to 1.83 (outer most). Curves 
are drawn at intervals of dX = 0.02. Outflow rate is maximum 
at some intermediate compression ratio. 
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these shocks is written. Here to compute sohd angles of the 
inflow and the outflow, we assume the half opening angle 
of the outflow to be 10°. Therefore, Qoui = 7r'^/162. 8i„ is 
given in the discussion following Eq. (3). In Paper I, the 
compression ratio R was assumed to be a parameter and 
no angular momentum was assumed a priori. Presently, 
we show the dependence on angular momentum. The gen- 
eral character, namely, that the outflow rate is negligible 
when the shock is weak {R ^ 1) and falls off gradually for 
strongest shock {R ^ 7), remains the same as in Paper I, 
however. There is a peak at about Rj^ ~ 2.8%. Note that 
for a given R, R^ increases monotonically with specific an- 
gular momentum A. This is because density of CENBOL 
rises with A. The curves corresponding to A = 1.71 and 
1.73 arc specially marked since there is a clear difference 
in tendency of the variation of Rm- For instance, below 
A ~ 1.72, very strong shocks are not possible at all and as 
a result the outflow rate has a lower limit. For A ^ 1.72 
such a limit does not exist. 

The general behaviour of the outflow rate can be un- 
derstood in the following way: when shocks are strong, 
they form very far out, and thus, even though the 
CENBOL area (which is basically the area of the base 
of the jet) increases, the net outflow rate is low. When 
the shock forms very close to the black hole, the temper- 
ature is high, and thus the outflow velocity is larger, but 
the CENBOL surface area goes down. Thus the product 
is low. For the intermediate cases the net effect is larger. 

For comparison with the analytical work presented in 
Fig. 2a, in Fig. 2b we present a similar diagram drawn 
using a numerical computation of the shock locations 
(Chakrabarti, 1989). Excellent agreement between these 
two figures implies that the approximations on which the 
analytical work was based are justified. All the features 
are reproduced well in Fig. 2a, except that for the weak- 
est shocks outflow rate is not as low as in the numerical 
calculation of Fig. 2b. 

We now present the nature of Rm when the outflow is 
also chosen to be adiabatic in Fig. 3. We used 6o/0i ~ 0.1 
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Fig. 2b: Same as Fig. 2a except that curves are drawn for the 
exact numerical solution. 
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Fig. 3: Ratio of the Outflow and the Inflow rates as a function 
of the compression ratio of the inflow when the outflow is 
adiabatic. The general nature of the function remains the 
same as that of the isothermal outflow. 

for reference. We observe that the peak is still located at 
around R =~ 4 and the outflow rate drops for very strong 
{R ~ 7) and very weak (i? ~ 1) shocks. We therefore 
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believe that our conclusion about the behaviour of Rm is 
generic. 



3.2. Two component advective flows 

Chakrabarti & Titarchuk (1995) proposed that the spec- 
tral properties are better understood if the disk so- 
lutions of sub-Keplerian flows are included along with 
the Keplerian flows. Recently, Smith, Heindl and Swank 
(2001), Smith et al. (2001), Miller et al. (2001) found con- 
clusive evidence of these two components in many of the 
black hole candidate accretion flows. While the matter 
with higher viscosity flows along the equatorial plane as 
a Keplerian disk (of rate Mk), sub-Keplerian halo matter 
(of rate Mh) with lower viscosity flanks the Keplerian disk 
above and below (Fig. 3a). Since the inner boundary con- 
dition on the horizon forces the flow to be sub-Keplerian, 
irrespective of their origin (Chakrabarti, 1990, 1996) mat- 
ter mixes (at say, r = rtr) from both the Keplerian and 
sub-Keplerian flows before entering a black hole to form 
a single component sub-Keplerian with an average energy 
and angular momentum of £ and A respectively. The spe- 
cific energy and angular momentum of the mixed flow is 
computed from: 

^ MkSk + Mh£h 



and 



A = 



MK + Mh 
Mk\k + MhXh 



(21) 



(22) 



Mk + Mh 

Here, £k, £h, and A/, are the specific energies and 
specific angular momentum of the Keplerian and the sub- 
Keplerian components at r = rtr respectively. 

Figure 4a shows a schematic diagram of the cross- 
section of a two-component accretion fiow. The transition 
radius (r = rtr) where the Keplerian disk becomes sub- 
Keplerian, and the shock location r = Ts, are indicated. 
Fig. 4b shows two solutions (marked I and II) of the equa- 
tions governing a two-component flow (Chakrabarti, 1996) 
where Xd/^K (Sub-Keplerian matter from the Keplerian 
disk) and Xh/Xk (Sub-Keplerian halo) are plotted as a 



Winds (Mout) 

r 




Sub ''^P'^f'sn 
Keplerian (^d'!^^) 




Fig. 4a: Schematic diagram of the cross section of two- 
component accretion flow. See text for details. 




1 2 

Radial Distance [log(r)] 



Fig. 4b: Solution of the two-component flow equations for 
two different viscosities. They axe merged to form a single 
solution as depicted in Fig. 4a. 

function of the logarithmic radial distance. Viscosities cho- 
sen for these two components are a = 0.04 and a = 0.01 
respectively. For r < rtr = 45 (lightly shaded region) 
the two sub-Keplerian flows mix to create a single com- 
ponent. For simplicity, we assume viscosity to be negli- 
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Compression Ratio (R) 



Fig. 5: Variation of outflow rates (left axis) with compression 

ratio at shocks (lower axis) . The upper axis gives the variation 
of sub-Kcplorian accretion rate and right axis gives the same 
for Keplerian accretion rate. 

gible in this region. Thus, the specific angular momen- 
tum and specific energy computed at r = rtr from Eqs. 
(21 & 22) remain constant (A) for r < rtr- Dark solid 
curve (marked III) shows the angular momentum distri- 
bution X/Xk of all possible mixtures of the two compo- 
nents which allow shock formation. We chose a case where 
Md + Mh = "^-OMEdd and vary the Keplerian component 
Md where Msdd is the Eddington accretion rate. 

In Fig. 5, the computed outflow rates are shown when 
the half opening angle of the outflow is 10°. In this case, 

~ \/ ^TT^ ^ , /„ . The left axis shows the rate of out- 

flow jfiout = Mout/ M^Edd as a function of the Keplerian 
disk rate (right panel) {md = Md/MEdd) and the halo 
rate (upper panel) {rhh = Mh/MEdd)- The lower axis 
gives the compression ratio at the shock. The most im- 
portant conclusion that can be drawn here is that the 
outflow rate steadily goes up as the Keplerian disk rate 
rhd decreases and the spectrum goes to a harder state. 
When the Keplerian rate is higher, the compression ratio 
is lower and the outflow rate is also lower. This conclusion, 
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drawn completely from dynamical considerations, is also 
found to be true from the spectral studies (CT95) where 
it was shown that the post-shock region cools down and 
the shock disappears (i? 1). Our work therefore hints 
that the outflow would be negligible in softer states. 



4. Discussion and Concluding Remarks 

CT95 pointed out that the centrifugal pressure-supported 
boundary layer (CENBOL) of a black hole accretion flow 
is responsible for the spectral properties of a black hole 
candidate. In this Paper, we present analytical results to 
show that this CENBOL is also responsible for the pro- 
duction of the outflows, and the outflow rate is strongly 
dependent on the inflow parameters, such as specific en- 
ergy and angular momentum. We showed that in general, 
the outflow rate is negligible when the shock is absent and 
very small when the shock is very strong. In intermediate 
strength, the outflow rate is the highest. As the specific 
angular momentum is increased, the outflow rate is also 
increased. This conclusion is valid when the flow is either 
isothermal or adiabatic. 

We also demonstrated how a realistic two-component 
fiow (TCAF) consisting of Keplerian and sub-Keplerian 
components produces a significant amount of outflow. 
Since matter close to a black hole is sub-Keplerian by 
nature, the two components must mix to form a single 
sub-Keplerian fiow which has positive specific energy and 
almost constant specific angular momentum. We showed 
that as the Keplerian rate of the disk is increased, the out- 
flow rate is decreased as the shock compression ratio ap- 
proaches unity. This conclusion, drawn from a dynamical 
point of view, is also corroborated by the spectral behavior 
as well — as the Keplerian rate is raised, the post-shock 
region is cooled due to inverse Comptonization and the 
shock disappears. This reduces the thermal pressure drive 
and the resulting outfiow rate is reduced. 
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